Aging in an infinite-range Hamiltonian system of coupled rotators 
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We analyze numerically the out-of-equilibrium relaxation dynamics of a long-range Hamiltonian 
system of A'^ fully coupled rotators. For a particular family of initial conditions, this system is known 
to enter a particular regime in which the dynamic behavior does not agree with thermodynamic 
predictions. Moreover, there is evidence that in the thermodynamic limit, when A'' ^ oo is taken 
prior to t —> oo, the system will never attain true equilibrium. By analyzing the scaling properties of 
the two-time autocorrelation function we find that, in that regime, a very complex dynamics unfolds, 
in which aging phenomena appear. The scaling law strongly suggests that the system behaves in a 
complex way, relaxing towards equilibrium through intricate trajectories. The present results are 
obtained for conservative dynamics, where there is no thermal bath in contact with the system. 
This is the first time that aging is observed in such Hamiltonian systems. 



At the very foundations of statistical mechanics, there 
are still some hypotheses whose validity rests merely on 
the extrapolation of observational facts and that have to 
be justified a posteriori. Among them, let us mention 
two assumptions that are intimately related to the issue 
addressed in this paper. The first one refers to the intro- 
duction of a probabilistic description of the evolution of 
a physical system. The second is related to the mechan- 
ical specifications that a system must fulfill so that the 
results of statistical mechanics can be applied [1] . These 
two points are closely related to the fundamental problem 
of establishing a connection between the dynamical be- 
havior of a system, described by the Hamiltonian H , and 
its thermodynamics. In that sense, statistical mechan- 
ics requires the existence of adequate conditions allow- 
ing to replace the dynamical temporal predictions by a 
probabilistic ensemble calculation that yields the correct 
equilibrium mean value of the relevant quantities. 

A very fast relaxation and a high degree of chaos and 
mixing are usually required in order to guarantee that a 
system orbit will cover most of its phase space in a short 
time. These questions have been extensively investigated 
for low dimensional systems, whereas, for extended sys- 
tems, with infinite degrees of freedom, the matter is still 
far from settled [2]. 

Due to the analytical and numerical efforts of many 
authors, it is today a well established fact that even 
very simple models, when analyzed in the thermody- 
namical limit iV — > oo, can yield results that bring to 
surface central questions about the foundations of statis- 
tical mechanics. In particular, in this work, we will con- 
centrate on an infinite-range (mean-field) Hamiltonian 
system that, despite its simplicity, exhibits a very pecu- 
liar dynamical behavior: depending on the initial prepa- 



ration of the system its evolution can get trapped into 
trajectories that will prevent the system from attaining 
equilibrium in finite time when the N oo limit is taken 
before the t ^ oo limit [3-8]. Unlike most models that 
display complex macroscopic behavior, this infinite-range 
model includes neither randomness nor frustration in its 
microscopic interactions. Furthermore, on one hand it 
can be exactly solved in the canonical ensemble, while 
on the other hand it can be efficiently integrated in the 
microcanonical ensemble. In that sense, it is an excellent 
starting point for analyzing the above mentioned basic 
questions. 

The system consists of N fully coupled rotators whose 
dynamics is described by the following Hamiltonian: 



) =K + V. (1) 



It is worth mentioning that this is a rescaled version of 
a nonextensive infinite-range Hamiltonian [9]. However, 
both of them share the same dynamical behavior after 
appropriate rescaling of the dynamic variables. 

In Fig. 1, we display the plot, T vs. U, where 
T = 2{K)/N is the temperature and U = {K + V)/N is 
the total energy per particle. The solid line corresponds 
to the canonical calculations [3] , which predict a second- 
order transition at Uc — 0.75. Above Uc constant specific 
heat is found and below Uc the system orders in a clus- 
tered phase. The symbols correspond to the numerical 
results obtained by integrating the equations of motion 
for N = 1000, 5000 rotators and until t = 1000. The sys- 
tem is initially prepared in a "water-bag" configuration, 
that is, all the angles are set to zero while the momenta 
are randomly chosen from an uniform distribution such 
that the system has total energy NU. By measuring the 
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nonequilibrium temperature (or equivalently the magne- 
tization) of a system started in these out-of-equilibrium 
initial conditions, one observes that, for a range of en- 
ergy values below the transition, the system enters in 
a quasistationary regime characterized by a mean kinetic 
energy that varies very slowly. Moreover, the value of this 
non-equilibrium temperature remains different from that 
predicted by canonical calculations. Actually, standard 
equilibrium is attained only after a time which grows 
with the size of the system, hence an infinite system will 
never reach true equilibrium [7,8]. In the quasistation- 
ary regime preceding equilibrium, trajectories are non- 
ergodic and the dynamics is weakly chaotic with Lya- 
punov exponent vanishing in the thermodynamic limit 
[5]. 

It is our objective here to show that the discrepancy 
between the results drawn from the dynamics and those 
derived from the canonical ensemble is closely associated 
to the presence of strong long-term memory effects and 
slow relaxation dynamics, a phenomenon usually named 
aging. Aging is one of the most striking features in the 
off-equilibrium dynamics of complex systems. It refers 
to the presence of strong memory effects spanning time 
lengths that in some cases exceed any available obser- 
vational time. Although aging has been seen in a wide 
variety of contexts and systems [10], some of them, ac- 
tually very simple ones [2,11], it is perhaps in the realm 
of spin glass dynamics where a systematic study of these 
phenomena has been carried out (see Ref., [10] and ref- 
erences therein). Systems that age can be classified into 
dynamical universality classes according to the sealing 
properties of their relaxation function. Moreover, these 
scaling properties contribute to a quantitative descrip- 
tion of complex phenomena, even in cases where a general 
theory is lacking [10,12]. 

Aging can be characterized by measuring the two-time 
autocorrelation function along the system trajectories. If 
the state of the system in phase space can be completely 
characterized giving a state vector x, then the two-time 
autocorrelation function is defined as follows: 



{x{t + tu,) ■ x{t^)) - {x{t + 1^)) ■ {x{tw)) 



(2) 



where at' are standard deviations and the symbol (• • •) 
stands for average over several realizations of the dynam- 
ics. In the case of a Hamiltonian system with N degrees 
of freedom, the state vector is decomposed in coordinates 
and their conjugate momenta, therefore we establish the 
following notation: x = {9,L). 

For systems that have attained "true" thermodynami- 
cal equilibrium, only time differences make physical sense 
when calculating relaxation quantities. In this case, it is 
expected that on an average the system will show only 
very short memory of past configurations. However, for 
systems exhibiting aging, a complex time dependence is 
observed in the behavior o{C{t + t.u,,tw), indicating long- 



term memory effects. In such a case, even at macroscopic 
time scales, the two-time autocorrelation function shows 
an explicit dependence on both times {t and tw) together 
with a slow relaxation regime. 

In order to integrate the motion equations numerically, 
we employed a fourth order sympleetic method [13] with 
a fixed time step selected so as to keep a constant value 
of the energy within a relative error AE/E of order 10~^. 
All the simulations were started from the water-bag ini- 
tial conditions explained above. 

In Fig. 2, we present the results of the numerical cal- 
culation of the two-time autocorrelation function (2) for 
U = 0.69. This value of the energy together with the 
water-bag initial conditions set the system into a par- 
ticular dynamical regime, in which ensemble discrepancy 
is more pronounced when finite size results are extrapo- 
lated to the thermodynamic limit. In the graph, features 
usual of aging phenomena can be distinguished. For a 
given tw the system first enters a quasiequilibrium stage, 
in the sense that temporal translational invariance holds, 
with C{t + tw,tw) « 1, up to a time of order tw After 
that, the system enters a second relaxation characterized 
by a slow power law decay and a strong dependence on 
both times. This phenomenology can be clearly seen in 
the curves for the largest i^'s. 

In Fig. 3, we show the best data collapse for the long- 
time behavior of the autocorrelation function, using the 
data of Fig. 2 corresponding to the three largest wait- 
ing times (U = 2048, 8192, and 32768). The resulting 
scaling law clearly indicates that for the whole range of 
values of t/tyj considered: 



C(t + t-u;, tw) — f(^~^ ) 



(3) 



where f{t/t^) ^ worth mentioning that 

this scaling is the same as observed experimentally in spin 
glass systems [14]. The values obtained for the scaling 
parameters are /3 « 0.90 and A « 0.74. In the inset, we 
exhibit an alternative representation of the data which 
yields a linear plot. It corresponds to \nq[C{t + tu],tw)] 
vs. t/t^, where the function lnq(a:), named gr-logarithm, 
is defined as follows [15]: 



lnq(a;) 



1-q 



(4) 



In this expression, q = 1 + 1/A, which for the data in 
Fig. 2, yields q « 2.35. Therefore, we can obtain a 
complete functional form of the autocorrelation function 
valid over the whole range of the scaling variable t/tf^ 
just by identifying the fimction f{x) in Eq. (3) with the 
inverse of lng(a;), that is, the g-exponential [15]: 



/(x)=e--= l-il-q)x 



(5) 



which naturally arises within the nonextensive statistics 
introduced by Tsallis inspired by the probabilistic de- 
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scription of multi- fractal geometries [16]. The same quali- 
tative behavior has been observed for other systems sizes, 
namely N = 500, 2000. 

This aging scenario contrasts with the time invariant 
behavior observed within the high energy phase, where 
no quasistationary regime is detected. In fact, let us dis- 
cuss Fig. 4 where we present the results of the calculation 
of the two-time autocorrelation function (2) for U = 5.0, 
well above the second order phase transition (i.e., in- 
side the homogeneous phase) , with water-bag initial con- 
ditions. What we observe here is essentially that the 
autocorrelation function depends on the two times only 
through their difference, that is, C{t + tyj,tyj) « C{t). 
Therefore, the presence of aging is related to the exis- 
tence of quasistationary states. It is important to em- 
phasize that, although the dynamics presents temporal 
translational invariance in the high energy regime, the 
relaxation of the system is very slow. 

There is nowadays growing evidence that aging is a 
very common dynamical phenomenon, associated to a 
great variety of physical systems. So far, there are two 
scenarios within which aging can emerge. On the one 
hand, the onset of aging in many systems derives from 
the presence of coarsening processes that give place to 
critical slowing down of the dynamics. In this case the 
scaling law of the two-time autocorrelation function is 
ruled by the following expression 

C{t + )^f{L{t)/L{t^)), (6) 

where L{t) is the mean linear size of the domains at time 
t [17]. On the other hand, aging also appears as a onse- 
quence of weak ergodicity breaking [10] and it is related 
to the complex fractal structure of the region of phase 
space that the system explores in time. This is the case, 
for instance, in the Sherrington-Kirkpatrick (SK) model 
and other spin glass models in which the complexity of 
the energy landscape is associated to a certain degree of 
randomness and/or frustration in the Hamiltonian. 

What is particularly remarkable in this work is the 
appearance of a complex aging behavior in a mean-field 
model lacking both randomness and frustration. Since 
the model wc arc analyzing in this paper is an infinite- 
range one, such as the SK model, a coarsening scenario 
has to be ruled out from the outset. Furthermore, all our 
results arc obtained in a conservative system without any 
thermal bath in contact with it. 

Moreover, it is worth noting that the scaling law found 
for the two-time autocorrelation functions below the 
transition, where a quasistationary regime is detected, 
points to a scenario very similar to that observed in 
spin glasses [14]. As occurs in spin glasses, there is 
weak breakdown of ergodicity which is consistent with 
the observation of weakly chaotic orbits, i.e., with a van- 
ishing Lyapunov exponent in the thermodynamic limit 
[5]. Drawing the analogy with spin glasses even further, 
our results seem to confirm that the system visits phase 
space confined inside very intricate trajectories (presum- 
ably nonergodic). This conjecture is also supported by 



features observed in fi space [6 8]. The fact that the re- 
laxation of the two-time autocorrelation function can be 
well fitted by a g-cxponential decay over the whole range 
of t/tyj deserves further investigation. Although a possi- 
ble connection with nonextensive statistics [16] is still not 
clear, wc believe that it would be interesting to examine 
this possibility. 

In summary, in this paper we have characterized the 
slow relaxation dynamics of a long-range Hamiltonian 
system through its aging dynamics. Our observation of 
the existence of aging in this Hamiltonian system and 
its characterization by scaling properties reminiscent of 
spin glasses is a result that can contribute to establish a 
unified frame for the discussion of the out-of-equilibrium 
dynamics of systems with many degrees of freedom. 
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FIG. 1. The full line corresponds to the canonical theo- 
retical caloric curve and the symbols correspond to numer- 
ical simulations for systems of size N = 1000 (circles) and 
N = 5000 (triangles), aX t = 1000, averaged over 50 realizar 
tions. Initial conditions are "water-bag" . 

FIG. 2. Two-time autocorrelation function C{t + tu,,tw) 
vs. t for systems of size N = 1000 and energy per particle 
U — 0.69. The data correspond to an average over 200 tra- 
jectories initialized in water bag configurations. The waiting 
times are t^=8, 32, 128, 512, 2048, 8192, and 32768. 

FIG. 3. Data collapse for the long-time behavior of the au- 
tocorrelation function C{t + tw,tm)- The data are the same 
shown in Fig. 2 for the three largest tw The gray solid line 
corresponds to eq{—0.2t/t'^^). Inset: luq-linear representa- 
tion of the same data, with q ~ 2.35. 



FIG. 4. Two-time autocorrelation function C{t+tw,tw) vs. 
t for systems of A'^ = 1000 and U — 5.0. The data correspond 
to an average over 10 trajectories initialized in water bag con- 
figurations. The waiting times are tw=8, 32, 128, 512, 2048, 
8192, and 32768. Inset: semi-log representation of the same 
data. 
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Figure 1 
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